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Direction-dependent Jeans instability in an anisotropic Bianchi type I space-time 
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We derive the metric for a Bianchi type I space-time with energy density that is dominated by 
that of a perfect fluid with equation of state p = wp and whose anisotropy is seeded by a fixed norm 
spacelike vector field. We solve for the evolution of perturbations about this space-time. In partic- 
ular, the Jeans instability in an expanding flat Friedmann-Robertson- Walker universe is modified 
by the presence of the vector field so that energy density perturbations develop direction-dependent 
growth. We also briefly consider observational limits on the vector field vacuum expectation value, 
m. We find that, if m is constant during recombination and thereafter, m < 10 14 GeV. 



I. INTRODUCTION 

The Universe is usually assumed to be homogeneous 
and isotropic on large scales. In this paper, we make 
progress in testing the assumption of isotropy by ex- 
ploring a particular model of isotropy-breaking involving 
fixed-norm spacelike vector fields minimally coupled to 
gravity. Asymmetries (like anisotropy) that have been 
very nearly ruled out in the current cosmological epoch 
could have been present in the early Universe; the sig- 
natures of asymmetries long ago could be very subtle 
and therefore not ruled out experimentally. Further- 
more, although rotational symmetry breaking in Stan- 
dard Model particle interactions is very tightly con- 
strained (see, for example, Refs. [H, H, H, 13, H, B 0])' 
isotropy-breaking through minimal coupling to gravity is 
not as well-constrained H, [T3, EH ■ 

In trying to understand a broken symmetry, it is often 
convenient to study a particular model that provides the 
mechanism for symmetry breaking. There is a class of 
natural models for breaking Lorentz invariance sponta- 
neously involving fixed-norm vector fields. Fixed-norm 
timelike vector field models that break Lorentz invari- 
ance but not rotational invariance have been extensively 

studied [H [H, M m EEES HE M M SH II M M 

[Hi HI]- For a review, see |27l |. 

Recently, there has been interest in the consequences 
of a small brea king of rotational invariance during the 
inflationary era [lllljp, M [H, HIH H|. A variant 
of the models of Refs. [I4jj7 |22j , and [12| provide a conve- 
nient framework with which to explore such consequences 
[1^, [36j]. The model involves a spacelike, fixed-norm vec- 
tor field, with Lagrange density [l4| . 
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Here A is a Lagrange multiplier that enforces the fixed 



norm constraint g^u^u" 



rn 



The Lagrange multi- 



plier term can be thought of as encoding the energy- 
minimizing effects of a smooth potential. The four-vector 
u a is spacelike and induces the spontaneous breaking of 
rotational invariance when m 2 > 0. (A timelike four- 
vector would break Lorentz invariance but not necessar- 
ily rotational invariance.) A model for the breaking of 
isotropy during the inflationary era with dynamics gov- 
erned by the action, 



S = 



d 4 x. 
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where p\ is a constant energy density, was considered in 
(29| in the interest of understanding effects of isotropy 
breaking during inflation on the cosmic microwave back- 
ground (CMB). 1 When coordinates are chosen such that 
the X3-axis is parallel to the direction along which the 
vector field gets its vacuum expectation value (VEV), 
the resultant space-time is given by [2^ |. 



ds 2 



-di 2 + a(<) 2 dxi + b(t) 2 dx 2 3 



where, 
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and, 
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(1 + 87rG/3im 2 )(3 + 32nGl3 1 m 2 ) 



(5) 



In the quest toward testing the standard assumption of 
spatial isotropy in our Universe, it is informative to con- 
sider the existence and properties of other anisotropic 
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1 A related 5-dimensional model was recently considered in [37 
a way to hide large extra dimensions. 
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space-times seeded by a fixed-norm vector field. To this 
end, we consider a universe with a fixed-norm spacelike 
vector field, permeated by a perfect fluid with equation 
of state p = wp where w is constant. We consider this 
perfect fluid case because (1) it can model relevant cos- 
mological epochs such as radiation (w = 4) and cold 
non-relativistic matter (w — 0) domination and (2) it is 
the next most mathematically tractable example of an 
anisotropic space-time. The case w = — 1 corresponds to 
a vacuum energy density as considered in [36| . 

We assume that the space-time geometry is approx- 
imately spatially flat; Ref. 38J found that an approxi- 
mately spatially flat geometry is preferred if anisotropy 
is to remain small. Considering this particular model 
of spontaneous isotropy-breaking when ordinary matter 
is present allows us to study the modified behavior of 
matter and metric perturbations — in particular the mod- 
ification of the Jeans instability. This particular model 
also provides an example of an expanding Bianchi type 
I space-time that does not evolve to be isotropic. Also, 
the deviation from isotropy can be tuned by a single pa- 
rameter. 

In Ref. [36[ , an analysis of vector field perturbations in 
flat space led to the requirement for energetic stability: 
Pi+fh+Pz = and f3\ > 0. This result carries over to our 
present analysis (since the analysis in [3(| was indepen- 
dent of gravity and other matter fields) . We shall further 
suppose that an analysis of energies in non-Minkowski 
space would lead to the requirement that the vector field 
Lagrange density take the form of a field strength tensor 
squared (£ u ~ F^ V F^ V ) in both Minkowski and non- 
Minkowski space-time. 2 Therefore, we consider, 

Cu = -\f^F^ u + X(u^ - m 2 ), (6) 

where F pv — V — V^u^. We have set (i\ — 1/2 in 
order to canonically normalize the kinetic term. Thus 
the parameter m 2 completely characterizes the deviation 
from isotropy. 

II. BACKGROUND SPACE-TIME 

We take the following action to govern our model, 
S = J d 4 x^j (j^jR + C fluid + CA , (7) 
where L u is given by Eq. ©. 



2 In Minkowski space, /3i + P2 + Pz = corresponds to a field 
strength squared Lagrange density, while in curved space a field 
strength squared Lagrange density corresponds to the more re- 
stricted case, fi\ = — /33 and j3% = 0. 



The energy-momentum tensor for w M derived from © 

is E3, 

T$ = F W F V " + g^C u + Wai (g va - ^f) V P F»° 

+ u u (g^ - ^f) V p FO° + ^ u ^ p FO°. (8) 
The fluid stress energy tensor is given by, 

Tjil l) = {p+piv^Vu+pg^, (9) 

where v p is the fluid's four-velocity and p and p are the 
energy density and pressure in the fluid's rest frame. Tak- 
ing the equation of state to be of the form p = wp, the 
stress-energy tensor takes the form, 

TtfP =p[(l+w)v„v l/ +wg fiV }. (10) 

We assume the fluid is homogeneous and isotropic in its 
rest frame, so p = p(t) . 

A. Solving Einstein's equations 

We choose coordinates such that the fluid is at rest. 
Since the fluid is isotropic, we still have the freedom to 
choose spatial coordinates such that the X3 axis aligns 
with the isotropy-breaking vector field . However, a pri- 
ori, the vector field may still have a timclikc component 
in the rest frame of the fluid. But, as is natural, we will 
assume that the vector field aligns such that it has no 
timelike component in the rest frame of the fluid. Wc 
hope to address this alignment issue in a future publica- 
tion. 

Taking the vector field to have no timelike component, 
we have, 

u M = m v / 333 7 7M3, (11) 

where <? M „ is the background metric. We make the ansatz 
for the background metric, 

ds 2 = -dt 2 + a 2 (t)dx 2 ± + 6 2 (<)(rfa; 3 ) 2 , (12) 

where rfx^ = (dx 1 ) 2 + (dx 2 ) 2 . This is an axisymmetric 
Bianchi type I geometry. Given Eqs. pip and (JSJ), the 
nonvanishing components of the background stress tensor 
for the vector field are [29| , 
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Note that the background solution satisfies the weak en- 
ergy condition if, 



There are two linear combinations of the field equa- 
tions (fl7|) - (fT9|) that are independent of p: 



l + w> 



4r 



3 + 4e' 



where we have defined, 



e = AirGm 2 . 



The field equations, 



G^ = 8ttG[T^+T^ 



(14) 



(15) 



(16) 



give the following non-trivial equations (p, = v = 0, p, = 
v = 1 and fj, = v = 3 respectively) for the background: 



a \ „ / a 
- +2 - 
a I \a 



8nGp(t) + 6 



- - \=-^Gwp{t)-e I 



iirGwp{t) 



(17) 



(18) 



(19) 



(2w + l) 



(23) 



2r 



2(2 

a 



(24) 



It is clear from the form of these linear combinations 
that the solution is of the form a(t) oc t a and b(t) cx t 13 . 
A third independent equation implies p(t) cx t~ 2 (when 
w ^ -1) and thus, by Eq. ([22]) . 



(2a + /3)(l+io) = 2. 
We find from these equations that, 

2(1 + 2e) 

OL — 

(l+tfl)(3 + 4e)' 
2 



= 



(l+tfl)(3 + 4e)' 



(25) 

(26) 
(27) 



The background metric satisfies the following equa- 
tions, 



Although it is not independent, it is helpful to consider 
the continuity equation, 



u— v I /i0 ' -'/iO J I background 



(20) 



We can simplify this equation in the usual way to deter- 
mine the relationship between the scale factors and the 
energy density p, 



-< 1+l °>is<» Jt >- < 21 » 



H a (*) = - =H b (t)(l + 2e), 
a 



Hb(t) -b-y (l + e)(3 + 4c)- 



(28) 
(29) 



These relationships between ff a , and p are iden- 
tical in form to the relationships, Eqs. (0, that were 
found in Ref. [29] for the case of a deSitter-like space- 
time (i.e. w = —1). For w ^ —1, we have just found 
that 



<>(h =,,(/„) (| 



(30) 



as in the case of a flat FRW universe. 3 However, since e > 
0, the direction in which the fixed-norm spacelike vector 
gets a VEV expands more slowly than in the isotropic 



The solution for the time dependence of p is given by, 

l+w 



p(t) _ fa 2 (t a )b(t a ) 



p(t ) V a*(t)b(t) 



(22) 



3 This indicates that evolution equations parameterized by cosmic 
time do not have a smooth limit w — ► — 1, since for w = —1 the 
energy density is constant. 
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case, while the transverse directions expand more quickly. 
We assume the deviation from isotropy will be small, so 
e << 1; in this limit. 



./ 8cl - ' V ' ' ' 



b ' 8(1 -!- ir)t V 1 3 £ 



(31) 
(32) 



In the isotropic limit (e — ► 0), the flat FRW metric for a 
perfect fluid with equation of state p = wp is recovered. 



III. CLASSICAL STABILITY 

A viable model of isotropy breaking must be at least 
classically stable against small perturbations. In [36j, it 
was shown that the space-time given by Eqs. ©-© is 
classically stable against small perturbations for a range 
of parameters. More specifically, the authors analyzed 
the equations of motion for vector field perturbations 
in flat space and showed that when the parameters in 
Eq. U) satisfy 0-y + 2 + 03 — and 0i > 0, the energy 
of field configurations that satisfy the Euler-Lagrange 
equations of motion is non-negative at all orders in per- 
turbations. When 0i + 2 + 03 ^ 0, negative energy 
modes propagate. Thus models with 0i + 2 + 03 
were deemed unstable. Also, Einstein's equations were 
solved in the limit where mode wavelengths are much 
smaller than the Hubble radius, and in the opposite limit 
where mode wavelengths are much longer than the Hub- 
ble radius. Modes in the short wavelength limit were 
shown to propagate as luminal plane waves with con- 
stant amplitudes. Short wavelength modes with superlu- 
minal phase velocities were found to propagate unless 
01+03 — and 02 = 0. Modes in the long wave- 
length limit were shown to have only constant parts 
and parts that decay as e -^H a +H„)t _ i/^{t)b(t)) or 

e- 167rGm2 ^ Hbt = b{t)/a{t). 

We now consider small perturbations about the back- 
ground given by Eqs. (33]), (©, (ESI), and J2S]). How- 
ever, let us first recall the story of perturbations in an 
expanding, flat FRW universe. For perfect fluids with 
constant equation of state, physical (i.e. gauge invariant) 
perturbations can behave differently depending on the 
scale of the wavelength of the perturbation. 4 Modes with 
wavelengths shorter than the Jeans wavelength (which is 
on the order of y/wH -1 ) will oscillate as plane waves. 
Modes with wavelengths much larger than H~ x decay or 



remain constant. Modes — in particular, the energy den- 
sity perturbation — at the intermediate scale can grow or 
decay; this is the Jeans instability in an expanding uni- 
verse. The instability is most efficient when w ~ since 
the Jeans wavelength goes to zero when w — > 0. 

We study only the very short wavelength (i.e. \y/wk\ 3> 
H) and long wavelength (i.e. \^/wk\ <C H) limits. How- 
ever, when w = 0, the very short wavelength limit as 
stated does not exist. For w — 0, we study the limit, 
k 3> H; the Jeans instability is manifest at this scale. 
Considering the very short wavelength and long wave- 
length limits is enough to establish stability as compared 
to the isotropic case. We assume that we can glean the in- 
teresting physics of the Jeans instability in an anisotropic 
FRW-like universe by considering an intermediate (short) 
wavelength scale in only the w — case. 



A. Parameterizing independent perturbations 

We define the ten (comoving) metric perturbations by 



ds 2 = -(1 + h Qa ) dt 2 + 2a(t) h 0t dtdx 1 
+ 2b(t) h 03 dtdx 3 + 2a(t) b(t) h^dx^x 3 



(33) 



+ a(t) 2 (Sij + hij^dxi + b(t) 2 (l + h 33 )(dx z ) 2 . 

We define the vector field perturbations and fluid velocity 
perturbations by, 



<5i/< = v 11 - ■ 



(34) 
(35) 



We make the approximation that the fluid's equation of 
state is unperturbed (so p—p — wp8 p ) and that the shear 
stress perturbations are negligible compared to the pres- 
sure perturbations. 5 This is an excellent approximation 
for a matter dominated era and a radiation dominated era 
before neutrino decoupling. (See, for example, (40l. l4lj|.) 

Since v^v^ = — 1, the timelike fluid velocity pertur- 
bation is given by, Sv° = — (l/2)/ioo, leaving only three 
independent fluid velocity field components. Also, the 
vector field constraint, u^u^ 1 = m 2 implies that Su 3 = 
— (m/26)/i33, so there are only three independent vector 
field perturbations. 

It is convenient to define the comoving fluid velocity 
perturbations (i.e. the peculiar velocity of the fluid) by 

(ie{i,2}) : 



Sv l 



V 1 



Sv J 



b(ty 



(36) 



4 Homogeneity of the space-time allows for perturbations to be 
Fourier-transformed in space and guarantees that modes with 5 Assuming Sp = wSp is equivalent to assuming that the pertur- 
different Fourier wave numbers will evolve independently. bations are adiabatic. See, for example Chapter 6 in Ref. [39l |. 
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Similarly we define the comoving vector field perturba- 
tions, 



Su° = mU°, Su 1 = m 



IP 

W)' 



(37) 



Finally, we define the fluid energy density contrast by, 



*„ = --!. (38) 



B. Gauge choice 

We can use diffeomorphism invariance to remove four 
of the seventeen perturbations just discussed. Under an 
infinitesimal coordinate transformation, — > x M + 
the vector field transforms to u x + Au x where, 

Au x = u»d^ x -ed ft u x , (39) 
and similarly — > + Av^ where, 

Av x = v"d^ x -ed„v x . (40) 
The metric transforms to g^ v + Ag^ v where, 

Ag^ = -(V m 6 + V^m)- ( 41 ) 

In analyzing short wavlength perturbations, we choose 
a gauge where the first order fluctuations in the (con- 
travariant) four-vector field about their background, 
= r) fl3 m/b(t), vanish (i.e. 8u^ L = 0). This gauge con- 
dition is satisfied by fixing, 



Au* 



1 st order 



u»d^ x - ed»u x 

171 ft it* _u m H 



(42) 



We shall refer to this as vector field gauge. Using this 
gauge allows us to carry over some results from (3|| since 
the vector field gauge was also used in that analysis. 

We find it more convenient to use a modified syn- 
chronous gauge, 



hoi 



and U 3 = 0, 



(43) 



in the analysis of long wavelength perturbations. This 
gauge condition is satisfied by fixing, 

-S(u 3 ) = Au x = -^d 3 e + e^H bl (44) 



b(t) 



b(ty 



-8{goi) = -a(t)h 0i = Ag ai 

= -dob - di£o + 2H a & 



(45) 



and, 



-<%03) = -b(t)h 03 = Ag 03 

= -dots - d 3 £o + 2H&. (46) 

The vector field satisfies the Lagrange multiplier equa- 
tion, u 2 — m 2 . This equation and either gauge choice 
imply that, 

(Sg^)u^ + 2g^(5u»)u v = h 33 m 2 = 0. (47) 



C. Perturbation equations 

By the twice-contracted Bianchi identity, conserva- 
tion of total energy-momentum density follows from Ein- 
stein's equations. Also, one can check that the equa- 
tions of motion for the vector field are equivalent to the 
conservation of the vector field stress-energy tensor. In 
general, Vj,TW 3 = identically at first order in pertur- 
bations, while the other three vector field stress-energy 
conservation equations are nontrivial. This — along with 
the twice-contracted Bianchi identity — means that there 
are only thirteen independent equations. The first order 
contracted Bianchi identities can be written as, 



(48) 



where, 



SE^ ee 8{G IIV ) - 8nG8(TM + T$). (49) 



See Appendix A for the full set of equations without 
any gauge condition applied. 



D. Short wavelength fluid dynamics and Jeans 
instability 

Einstein's equations in the very short wavelength limit 
(\/wk ^> aH) are independent of p and the background 
metric. Thus, since we use the same gauge as in [36| . 
the solutions for metric perturbations in the short wave- 
length limit found there (and described at the beginning 
of this section) apply also to the present case of a perfect 
fluid. The fluid pertubations must have the same plane- 
wave form as the metric perturbations when the metric 
perturbations are non-zero. The fluid stress-energy con- 
servation equations relate the metric perturbation ampli- 
tudes to the fluid perturbation amplitudes. 

In principle, there can also be a non-gravitational fluid 
perturbation mode governed by the fluid equations in 
flat space. The first order perturbed fluid stress-energy 
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conservation equations (in flat space) are: 6 

0= (l + u;)V- V + dtS p , (50) 
= wV5 p + {l + w)d t V. (51) 

We may decompose V as: 

V = V$ + £> ; V-B = 0. 

Then taking the time derivative of Eq. (|50|) and the 
divergence of Eq. (|5"Tj) and substituting in for V • V gives, 

(raV 2 - d 2 t )5 p = 0, (52) 

while taking the gradient of Eq. (|50|) and the time deriva- 
tive of Eq. (f5Tj) and similarly substituting gives, 

(wV 2 - <9 t 2 )V$ - d 2 B = 0. (53) 
The solutions for $ and S p are plane waves, 



d p = A p e luJt±lk - s and V$ = A$e ia,i±ife ' £ , (54) 



where, 



oj 2 = w(fc ■ fc) and A$ 



w k 



1 



|fc| 



(55) 



The divergence-free part of the fluid velocity, B, must be 
independent of time. We recover the standard result: the 
speed of sound of the perfect fluid is y/w. 

In the case of a cosmological constant (w = —1) the 
energy density perturbation is independent of space-time 
coordinates and thus remains a constant. 

Special attention must be paid to the case w = 0. It 
was assumed that ad t ~ k in the short wavelength analy- 
sis of [36| , but Eqs. (|54|) and (j55|) show that this is not a 
valid assumption for modes with non-zero energy density 
contrast. If <9 t ~ H ab , then none of the ten Einstein's 
equations, 8E pvl nor the fluid energy-density conserva- 
tion equations (see Appendix A) constrain 5 P in the limit 
k H when w = 0. To find the evolution of 5 pi we must 
find a fc-independent linear combination of the equations 
that involves S p . We have found such equations in the 
cases *3 = and k± = 0. 

When k± ^ and k% = 0, one finds the following 
fc-independent equation, 



--(3 + 4e)Jtf + 2H b (l + 2e)dt + df ) 8 P 



~(3 + 4e)ff b 2 - 1(3 + 4e)i7 b a t j /,„„ = 0. Cli! 



For a review of hydrodynamics in the absence of gravity see, for 
example, [431 • 



In the short wavelength limit, each of equations (|A6[) and 
(|A10[) tells us that /iqo vanishes; thus for this mode, 



S p + 2H a S p 



4ttG 



P Sp. 



(57) 



,(l + e) 

The solution of this equation is, 

s p = c 1 t 2 ^ 3 +^ +c 2 t-\ 

Note that t 2 /(3+4<0 « an d t -i ^ 

When = and ^ ^ 0, we find the ^-independent 
equation, 

(-(3 + 4e)(l + e)H 2 + AH b d t + 2d 2 ) 5 P 

+ (-(3 + lOe + 4e 2 ) J ff h 2 + 2(1 + e)H b d t ) h 00 

1 
2 



3etf 2 + (3/2 - e)H b 8 t + ^-d 2 ) h u = 0. (58) 



In the short wavelength limit when k± — 0, Eq. (|A11[) 
implies that /loo vanishes while Eq. (|A3|) implies that hu 
vanishes. Therefore, for this mode the density contrast 
evolves as, 



S p + 2H b S p = (4nGp)S p . 
The solution of this equation is, 

6 p = Citi++C 2 ti-, 

where, 

-1 + 4e± 5Vl + 3e(l + e)(4/5) 2 



(59) 



7± 



2(3 + 4e) 



1±5 
~~6 
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-4±- , 



In the isotropic limit, e — > 0, 



2HS P = (4nGp)5 p , 



(60) 
(61) 

(62) 



where H = lim e ^o H a b and thus the Jeans instability for 
a flat FRW space-time is recovered: S p ~ Cit 2 / 3 + C 2 t~ 1 . 
The above solutions tell us that, in the anisotropic case, 
density perturbations in directions transverse to the di- 
rection in which the vector field gets a VEV grow more 
slowly while those modes with wavenumber in the par- 
allel direction grow more quickly than in the isotropic 
case. 



E. Long wavlength limit 

Having completed an analysis of short wavelength per- 
turbations, we now consider the long wavelength limit of 
the perturbation equations. 
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Recall that H b is a function of time. In the modified 
synchronous gauge, we have the following field equations 
in the long wavelength limit, 



SEo 
1 + e 



= -(3 + 4e)H 2 (S p + /loo) + H b d t h 00 = 0, (63) 



a(t) 



8E, 



03 



b(t) 



= {I + w){'i + Ae){\ + e)HlV l = 0, (64) 



= (l + w)(3 + 4e)(l + e)H b 2 V 3 = 0, (65) 



8E, 



2 = <%i?b(l + e)[-w(3 + 4e)H b {8 p + h 00 ) + d t h 00 ] 



o(t) 

+ 2 e '(i e J)t (( 3 + 4e )#^ + dt) hk = 0, (66) 



28 E. 



(3 



= e(3 - u;(3 + 4e))H 2 [(l + 4e)/i i3 + (2 + 4e)tT] 



a(t)6(t) 

+ (1 + 4e) [(3 + Ae)H b d t + d 2 ] h a 
+ 4e[4(l + e)H b d t + 2 }U % = 0, 



(67) 



33 



b{tf 



= H b (l + e)[-w(3 + 4e)H b {8 p + h ao ) + d t h 



DO 



-((3 + 4e)ff b dt + d 2 )/^ = 0. 



(68) 



The fluid equations of motion in the long wavelength 
limit are, 



8(V^ ) 
P 

(1 + w)pa(t) 
8(V^ 3 ) 



dt8 P + -(w + l)d t hu = 0, 



(69) 



[(1 + 2e - w(3 + 4e))H b + d t ]V l = 0, 

(70) 

[(l-w(3 + 4e))H b + d t ]V 3 = 0. (71) 



(l + w)pb(t) 

And the equations of motion for the vector field are, 



8 {P „V P FP°) 



= m 



1- 3w 



+ 2e(l-w) 



U° = 0, 
(72) 



S(P i(T V p FP° 



ma(t) 



= [H%ew(3 + 4e) - d 4 2 ] {U l + h 



;3 J 



ei/ h 2 ([/ I -/i i3 )-(3 + 4e)^a t ?7 I 



(2 + 46)^^=0, (73) 



where = 3 MCT 



The coupled fields fall into two uncoupled sets 
{8 P , hoo , ha} and {C/ 1 ,^}, while the remaining degrees 
of freedom are uncoupled. The solutions take the form, 



8 P (t) = -c p -c 00 1 -±^(aHY-^, 



h 00 (t) = c p + c oa (a 2 b)( w - i y 2 



hu(t) = h 2 2(t) + C 2 2 = - 



1 



1 + w 



8 P (i) + en 



h 12 {t) = c 12 + d 12 (a 2 b) ' , 
M*) = J(cS ) +«S ) (aV 



U 



where, 



4e^ 



(74) 
(75) 
(76) 
(77) 

(78) 

-1 

(79) 
(80) 

(4) (u>-l)(3 + 4e)(l + 4e) 



sW 1 



843 ) = 4 2) (6-(3 + 4 e )( w + l)), 



((3 + 4e)(u>- 1) + 2) 



The constants {c p , coo, •••} are fixed by initial conditions. 
The decaying parts of all comoving perturbations except 
his and U l go as (a 2 ^)"^. Special attention must be 
paid to ha and U % because each has a term proportional 

tO a/b = f4e/(3+4e)(l+™) which growg w[th t j me 

Consider a gauge transformation of the form (x^ — > 



£ M ) where, 
Co =6=0 



and 



&(t, 2) = c t za(ty 



(82) 



By Eqs. (|3"9")) - (|4"Tj) . this gauge transformation leaves all 
perturbations except and U l unchanged. In partic- 
ular the modified synchronous gauge conditions remain 
undisturbed. The gauge transformation simultaneously 
changes U l and ha by, 



Ah 



(3 



a ■ a 

-c,- and At/ = Ci T . 




(83) 



Thus we see that the growing part of both U % and ha 
can be gauged away because of the relative sign of the 
terms multiplying a/b in (|78p and (|79p . 

Therefore, for w < 1 (true for ordinary matter) all 
modes with wavelengths larger than the Hubble radius 
decay to a constant. In particular, h^ and U l decay to 
zero on a characteristic time scale proportional to 1/e. 

We recover the generic behavior of flat FRW pertur- 
bations in the long wavelength limit. As stated at the 
beginning of the section, FRW perturbations on scales 
much larger than the horizon have constant and decay- 
ing parts. 
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For good measure, let us consider the exact correspon- 
dence in the w = case. For nonrelativistic, pressureless 
(w = 0) matter, the metric is given by (|12p with, 



a = b oc 77 oc t 



2/3 



(84) 



where rj = Jdt/a is the conformal time. One can pa- 
rameterize scalar perturbations to the FRW metric by 



ds 2 = a{rif 



- (1 + 24>)dif - 2B tin dx m dr]+ 
{{l-2^)8 mn -2E >mn )dx m dx r ' 



where to, n S {1, 2, 3}. Then the function, 



1 



$ = ~(a(B-E')Y 
a 



(85) 



(86) 



where ' denotes -g-, is gauge invariant. The relationship 
between the comoving perturbations, /i^„, defined in (|33[) 
and the perturbations defined in (|85p in Fourier space is, 



^ 2 



h h 



s, 



E 



1 



k 2 



run I hrnn : 



'■'ill n \ h mn , 



(87) 



B = i^-ho r , 
k 2 



hoo 



(89) 



In our modified synchronous gauge, B — 0. So, in the 
isotropic limit (b — > a) when w = 0, we find the Fourier 
transformed gauge invariant function to be, 

-j- h 00 1 f k m k n x \ ,a' , „ 

* = —fZ I d m n ( — ft mn + h mn ), 

2 4fc 2 V k 2 J a 



coo J?o 



1 Vo ( ^k m k n 
~2 ^ ~2~^ ik 2 ^ \ k 2 
x (2loC(»)+')oCW), 



rf' 



£oo 
2 



1^1 2 f r^k m k n 



x (2']o'»™(')o)+1oC('?o)) 



(90) 



since a oc ?7 2 . Here we've defined 0(770) = 1. In the long 
wavelength limit, \kr]\ << 1, and so the second term in 
the brackets dominates over the first for all times. Thus 
we have, in the long wavelength limit, 



$ ~ d(x) + C 2 (x)V~ 



(91) 



This is precisely the functional form found for $ in a 
w = flat FRW universe (39|. One can work out similar 
correspondences for other gauge invariant quantities in 
the isotropic, long wavelength limit. 



IV. OBSERVATIONAL LIMITS 

The isotropy of the microwave background was used 
to place present-day limits on the density of matter with 
general anisotropic stresses in (28[. 7 . If we assume that 
the VEV of the vector field, to, was constant from the 
time of recombination, then, following [28| . we may cal- 
culate the temperature anisotropy due to the direction- 
dependent evolution of the photon temperature since re- 
combination, 



T Xl = T X2 = Tq—— = T exp 
a(t) 

T X3 = T o^y = T oexp 



H a dt 



- / H h dt 



by, 



ST _ T X1 - T X3 



T 



T 



exp 







-(l + 2e) / H b dt 
2< ( - ^ H b dtj exp 



o «o , ao 

2t W) g \W) 



(92) 
(93) 

(94) 



exp 



H h dt 



H h dt 



Thus, 



ST, 



y I today 

Observations show that [II] , 



« 2e(l + z rec ) log(l + z rec ). 



% <io-. 



(95) 



(96) 



(97) 



Taking 1 + z r 
AnGm 2 , 



1100 then we have the limit on e 



2e < 10" 9 . (98) 

Plugging in (%-kG)- 1 / 2 = 2.43 x 10 18 GeV gives the limit, 

to < l0 14 GeV. (99) 

Presence of the vector field during inflation also af- 
fects the CMB. The analysis in [Hj], indicates that the 
parameter defined in [29l | , g* , characterizing deviation of 



7 An analysis of the of the cosmological evolution of matter sources 
with small anisotropic pressures in various cosmological epochs 
was carried out in [28l An anisotropic equation of state for dark 
energy was considered in [43l |. 
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the power spectrum from isotropy during the inflationary 

era, 



P(k)=P(|k|)(l + 5 *(k-n) : 



(100) 



can be measured by the Planck experiment down to a 
value of (roughly) 10 -2 . The parameter is roughly 
equal to e = AirGm 2 [2{|. Thus, the vector field could 
have had a much larger VEV during inflation than during 
recombination (and thereafter) and still escape detection 
by Planck. If, however, the vector field has a constant 
VEV throughout cosmological history, the limit (|9"9"|) in- 
dicates that the effect of the vector field during inflation 
will be undetectable. 

Bounds from primordial nucleosynthesis were also con- 
sidered in 28]. Primoridial nucleosynthesis limits the 
change to the mean expansion rate imposed at the epoch 
of neutron-proton freeze-out. If the rate of cosmic expan- 
sion during Big Bang nucleosynthesis is slower, then weak 
interactions freeze out later, at a lower temperature. A 
lower temperature at freeze-out results in a decrease in 
the primoridal 4 He to H mass ratio, which is constrained 
by observation. It was noted in [28| that if the anisotropy 
does not decay (or decays very slowly) between the period 
of freeze-out and matter-radiation equality, then bounds 
from the temperature anisotropy are actually stronger 
than bounds from primoridal nucleosynthesis. For ex- 
ample, a bound from primordial nucleosynthesis (due to 
the difference between the effective Newton's constant 
measured in the solar system versus the effective cosmic 
expansion Newton's constant) on the VEV of a similar 
timelike vector field was found; for a normalization sim- 
ilar to ours, found m < 10 18 GeV. 

It would also be interesting to consider how a fixed- 
norm spacelike vector field modifies gravity in the weak 
field limit. Furthermore, it would be interesting to see 
whether the direction-dependent Jeans instability that 
we found significantly affects structure formation. 



where e = AirGm 2 . 

We went on to consider first order perturbations about 
the background space-time, assuming negligible shear- 
stress, and found that the space-time is as perturba- 
tively stable as the flat FRW universe. We derived an 
anisotropic generalization of the Jeans instability for an 
anisotropic w = FRW-like universe. We showed that 
energy density perturbations parallel to the direction of 
the vector field VEV are governed by Eq. (|5T 



2H b S p = {4nGp)5 p , 



(104) 



while the perturbations perpendicular to the VEV are 
governed by Eq. (|59|) , 



5 p + 2H a 8 p = {A^Gp/{\ + e))5 p 



(105) 



These equations imply that the (parallel/perpendicular) 
perturbations grow at a (faster/slower) rate than the 
usual Jeans instability predicts for a flat FRW universe. 

We computed the rough limit, m < 10 GeV, assum- 
ing the vector field was present during recombination and 
thereafter. The effects of the vector field, if present only 
during inflation, are much harder to detect. Stronger 
bounds might arise from analyzing solar system tests of 
gravity, structure formation, or SNla redshift data. 
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APPENDIX A: PERTURBATION EQUATIONS 

We Fourier transform the metric perturbations in 
space, 



V. CONCLUDING REMARKS 

In this paper, we considered an anisotropic Bianchi 
type I space-time with energy density that is dominated 
by that of a perfect fluid with equation of state p = wp 
and whose anisotropy is due to a fixed-norm spacelike 
vector field. We derived the background space-time 
(Eqs. (HU), (Ell) and (2^), 



ds 2 = -At 2 + a(i) 2 dx 2 + b(t) 2 dx 2 3 , 



H a (t) 



H b (t) = I 



H b (t)(l + 2e), 

I 8wGp{t) 
V(l + £ )(3 + 4e) : 



(101) 
(102) 

(103) 



h,j, v {t,k) = / d 3 xh pl/ (t,x)e 



(Al) 



as well as the energy density contrast, vector field and 
fluid velocity perturbations. Because the background 
space-time is spatially homogenous, perturbations with 
different comoving wavenumbers will evolve indepen- 
dently. The physical wavenumbers £^3 are defined in 
terms of the comoving wavenumbers by, 



and 



_fc3_ 

b(ty 



(A2) 



The following are the Fourier transformed first-order 
field equations, SE^ = SG^ — SttGST^ = 0, for 
the metric, vector, and fluid perturbations defined 
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in Eas.(|53 > ([3"5 ]) . We use the background solution of {1,2} and k 2 = k\ + k 2 . 
Eq. (|28|) and Eq. ([29]) to eliminate H a and p. In these 
equations, repeated {i,j, k,l} indices are summed over 



SE 00 = -(3 + 4e)(l + e)H^(S p + h 00 ) - 2ieH b k 3 U° - 2(1 + e)H b (ikjh oj + ik 3 h 03 ) + \{S l3 k 2 - k, I. ,)''., 

- kjhhja + (1 + e)H b d t h n + -{kl + 2(1 + e)H b d t )h 33 = 0, (A3) 

SE 1 ~* 

' •" !H = (1 + w)(3 + 4e)(l + e)H 2 (V l + h oi ) + 2ieH b k 3 U l + (1 + e)H b ihh 00 + -(S l3 k 2 - hk^Oj 



a(r) 2 

- 7^^3/103 + ~«&3 (2ei?b + <9 t ) h l3 + ^(ikjdthtj - ikidthjj) - ^ik~id t h 33 = 0, (A4) 



— ^ = (1 + w)(3 + 4e)(l + e)H 2 (V 3 + h 03 ) - 2ek 2 U° + 2e(l + 2e)^LT + 2eik l d t U l + (1 + e)H b ik 3 h 00 
b[t) 

- ^kik 3 h 0l + ^((1 + 4e)P - fc|)/i 03 + ~(1 + Ae)ik {2eH b + d t )h j3 - -ik 3 d t hu + eik 3 d t h 33 = 0, (A5) 



-«7(1 + e)(3 + 4e)H b 2 S l3 (5 p + h 00 ) - 2ieS l] H b k 3 U° + ~ (kikj - <5y (P - 2(1 + e)H b d t )) h oa 



a(t) 2 

+ ik^ (2(1 + 2e)H b + d t ) h 03 - e l{l e j)k (ik t (2(1 + e)H b + d t ) h ok + k 3 kih k3 ) 

+ (kl + (3 + ie)H b d t + Of) h lk + hfakj - 8 l3 kl)h 33 - <^((3 + Ae)H b d t + d 2 )h 33 = 0, (A6) 

25E 

= -ik 3 (2H b (l + e) + d t )h 0l - (1 + Ae)iki{2H b + d t )h 03 - k 3 (kjh i: j - hhjj) + hk 3 h o - k 3 h i3 - 2ehk 3 h 33 

+ (1 + 4e)(P<5y - kikj)h j3 + (1 + 4e) [(3 - w(3 + 4e))ei? b 2 + (3 + 4e)H b d t + d 2 ] h i3 + 4eik z (2H b + d t )U° 
+ 2e ((1 + 2e)(3 - w(3 + 4e))# h 2 + 8(1 + e)H b d t + 2d 2 ) U l + 4e(%P - hk^W = 0, (A7) 

(JJJoo 1_ 1 

— _ = - w {3 + 4e)(l + e)H^(S p + h 00 ) - -k 2 ± h 00 + -{kikj - k^^h.j - 2ehk 3 h i3 + ifc;(2(l + e)H b + d t )h oi 

- 2ik 3 e(2(l + 2e)H b + d t )h 03 + (1 + e)H b d t h 00 - ^((3 + Ae)H b d t + d 2 )h u + 4e(l + 2e)H b ik 3 U° - 2ehk 3 U l 
+ ek 2 ± h 33 + 2iek 3 d t U° + e((3 + 4e)i/ b 9 t + a t 2 )^ 3 3 = 0. (A8) 
We also have the following equations resulting from the fluid equations of motion, 

^ = -(1 + w)ihV l - (1 + w)ik 3 V 3 - d t S p --(w+ l)d t (h u + h 33 ) = 0, (A9) 

P " 
S(X7 tC/Om. 1 ) _ 1 
^ 11 ^ l) = wikiSp - (1 + w)((3w - 1 + 2(2w - l)e)H b - d t )(V % + h oi ) + -(1 + ^)ifc i /i 00 , (A10) 

V ^ ^ = ^fcs^p - (1 + w)((3w - 1 + 4we)H b - d t )(V 3 + h 03 ) + -(1 + w)ik 3 h m . (All) 
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The equations of motion for the vector field are (P^a = Q\ux — "^a" ) : 

HPoa Z pFn = - h " Q( j - 3w ) + ^ - + (*i + ~ k ^ u ° - ^ ( ai2 ) 

- + 2e)tf b + a f )(C/ 1 + + -H b ik 3 (h Q0 - ha) - -ik 3 d t h S3 + H b ik t h l3 = 0, (A13) 

1 - & L = -ih{H b + dt)(U° - M + + h jS ) - ik 3 H b h 0l - {2eH 2 b - H b d t )h l3 

ma(t) 

+ \hhh 33 + {Htell + w(3 + 4e)) - k 2 ± - kj - (3 + Ae)H b d t - df)(W + h i3 ) = 0. (AH) 

I 
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